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We develop the Lagrangian perturbation theory in the general relativis- 

'6<\ 

0> tic cosmology, which enables us to take into account the vortical effect of 

a^ . 

' the dust matter. Under the Lagrangian representation of the fluid flow, the 
^ '. 

T-H I Based on this, the coupling between the density and vorticity is clarified in 

I a non-perturbative way. The relativistic correspondence to the Lagrangian 

in ; 

' perturbation theory in the Newtonian cosmology is also emphasized. 

o ■ 
^ ■ 
o : 

■ I. INTRODUCTION 

The dynamics of fluids in the expanding universe is of great importance in the cosmol- 

o ■ 

-t— > ' ogy. For investigating such a dynamics, the Newtonian treatment is often used as a good 
■ 

' approximation for the region l/L <^ 1, where / is the scale of fluctuations of fluids and L 



^ ■ corresponds to the Hubble radius IQ. In the Newtonian cosmology, the Lagrangian pertur- 
bation theory in the dust universe has been developed for non-linear density fluctuations, up 



to the caustic formation In order to take into account the relativistic correction to the 

Newtonian dynamics, the cosmological post-Newtonian approximation has been formulated 
Furthermore, in the cosmological post-Newtonian approximation, the Lagrangian per- 
turbation theory has been also discussed Such a treatment beyond the Newtonian 
approximation may become important because of not only the theoretical interest but also 
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the recent progress in the observational cosmology. For example, the Sloan Digital Sky Sur- 
vey (SDSS) aims at the deep survey over several hundred Mpc |[T0|, in which it is not clear 
whether the Newtonian treatment is sufficient for such a large area. 

Relativistic theories of linear perturbations have been developed |l],|Tl],|T2| . The second- 
order extension of relativistic perturbation theory has been also developed |T^. However, 
they still depend on the assumption that the density fluctuation is small. In order to over- 
come the drawback, several perturbative approaches for the nonlinear dynamics have been 
proposed, which are based on the fluid flow approach [l^, the gradient expansion method 
p!5Hl7|, or the relativistic Lagrangian approach It is shown that the relativistic post- 

Zel'dovich approximation, a second-order extension of the relativistic Lagrangian approach, 
successfully express the non-linear evolution of the density contrast with higher accuracy 
than the conventional second-order theory ||2CI|| . 

So far, however, these relativistic approaches are all restricted within the irrotational 
dust. There has been no established way of calculating non-linear evolution of the density 
fluctuation without such limitation. Our paper is aimed to add to knowledge in literature 
by presenting a relativistic framework of such calculation. 

This paper is organized as follows. In section II, we show the integrals of the density 
and the vorticity in the general relativity, in the simple manner using the Lagrangian con- 
dition. This fact suggests strongly that the Lagrangian condition allows us to formulate 
the Lagrangian description in the relativistic cosmology. Under this condition, section III 
presents a perturbative Lagrangian approach. Summary and Discussion are given in section 
VI. For comparison, the vorticity in the Newtonian cosmology is discussed in the appendix 
A. Residual gauge freedoms in the Lagrangian condition are clarified in the appendix B. 
Greek indices run from to 3, and Latin indices from 1 to 3. We use the unit, c = 1. 
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II. LAGRANGIAN DESCRIPTION OF DUST FLUID 



Let us consider a simple dust universe, in which the energy momentum tensor is written 

as 

T^"^ = pu^'u". (2.1) 

The conservation law T'^'^^^ = gives 

P;X + P<;. = 0, (2.2) 
u'^.^u" = 0, (2.3) 

which are called as the continuity equation and the geodesic equation, respectively. 
The vorticity uj'^ of the fluid flow is defined by [pT| 



^ l^a/..p^^^^.^^ (2.4) 

where e"^'^/' denotes the complete anti-symmetric tensor with e"^^^ = 1/ \/—g and g = 
det{g^tj). From the geodesic equation Eq. ( |2.3| ), we obtain the propagation equation for the 



vorticity 21 



cu'^j^M^ + m'^.^cj^ = u^.^to". (2.5) 

Using Eq. (|2.2|) , we have 



(^),.."^.%(^), (2.6) 

which may be called as the relativistic Beltrami equation (cf. Eq. ( |A10|) in Appendix A). 
The Einstein equations are decomposed with respect to the fluid flow as follows: 

G^^uV = SnGp, (2.7) 

G^^u^P--^ = 0, (2.8) 

G.^P^^P'^f, = 0, (2.9) 

where P'^ is the projection tensor 



(2.10) 



So far the treatment is fully covariant. In the following, we adopt the Lagrangian condi- 



tion (e.g. ||22|), in which the components of the matter 4- velocity take the values of 



m'^ = (1,0,0,0). 



(2.11) 



Under this condition, we immediately have goo = —1 and = {—l,goi)- Furthermore, 
since the geodesic equation Eq. (|2.3|) under the Lagrangian condition simply tells m^_o = 0, 
Ui{= goi) are functions of spatial coordinates only: 



Ui = Ui{x), and goi = goi{x). 



In the Lagrangian description, the continuity equation (|2.2|) is simply 



(Pv^),o = 0. 



(2.12) 



(2.13) 



Therefore, 



p{x,t) 



9{x,to 



-p{x,to 



(2.14) 



Under the Lagrangian condition, the determinant of the metric tensor can be expressed as 



5, = -(1 + Y^goigoj) det{gij), 



(2.15) 



where 7*-^ is the inverse of the spatial metric gij. The relativistic Beltrami equation ( p.6| 
also becomes simply 



which is integrated to give 



— ,0 = 0, 
P 



P~ P to 



(2.16) 



(2.17) 



This is also expressed as 
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uj''{x, t) 



^^'''^'^u;\x,to). (2.18) 



The component is not independent of cu*. Using the relation cu'^m^ = 0, we obtain 



a;° = (7o.c^\ (2.19) 



The result Eq. ( p.l7|) tells us that the vorticity is coupled to the density enhancement 



and vice versa. In particular, if the vorticity does not vanish exactly at an initial time, 
the vorticity will blow up as the density grows larger and larger (i.e. in the collapsing 
region), even if it has only the decaying mode in the linear perturbation theory. It should 
also be emphasized that our results Eqs. ( 2.14 ) and ( p.l7|) in the fully general relativistic 



treatment precisely correspond to those in the Newtonian case (see Eqs. ( |A7| ) and ( [A13| ) in 
the Appendix). 

III. PERTURBATIVE LAGRANGIAN APPROACH 

In the previous section, we solved exactly the equations for the density and the vorticity. 
The results Eqs. ( |2.14 ) and ( p. IS ) show that p and u;* are completely written in terms of 



the determinant of the metric tensor and their initial values. In this section, we solve the 
metric perturbatively. We assume that the background is spatially flat Friedmann-Lemaitre- 
Robertson- Walker (FLRW) universe. The extension to the spatially non-fiat case must be 
a straightforward task. The perturbed metric is decomposed into 

goi = B^i{x) + bi{x), (3.1) 
gij = {6ij + 2HL5ij + 2HT,ij + (/ijj + hj^j) + 2Hij) , 

where B, Hi, and Ht are scalar mode quantities, bi and hi are the vector (transverse) mode, 
and Hij is the tensor (transverse-traceless) mode satisfying 

b\i = 0, (3.2) 
h\i = 0, (3.3) 



H\ = 0, 
= 0. 



(3.4) 
(3.5) 



Raising and lowering indices of the perturbed quantities are done by 6'^^ and 6ij. Using the 
perturbed metric Eq. (|3.1| ), the general expression for the perturbed Einstein tensor up to 
the linear order is as follows: 



(3.6) 



ct a a 



(3.7) 



5,: 



It ~\~ S—Hx 
a 



+ 



[h, + 3-hi - —bAj + ihj + 3-hj - —bj 



(3.8) 



where an overdot (') denotes d/dt^ = S^^didj, and we have used the fact a{t) = t^^^ so 



that G^uP iP^j does not have the background quantity. 



A. Scalar perturbations 

Using Eqs. (|3.7| ) and (|3.8|) , the Einstein equations for the scalar perturbations are 

{HL + ^By=0, (3.9) 

Hl + 3-Hl = 0, (3.10) 

a 

Ht + S-Ht = \ (Hl + -B). (3.11) 

In order to solve the above equations, we use the residual gauge freedom (see Eq. ( [B^ in 
Appendix B) to set 5 = 0. Then, from Eq. ( |3.9| ) we have 

Ei^ = Hl{x) = ^^^{x). (3.12) 



Next, from Eq. (|3.11|) , we obtain two independent solutions for Ht'- one is proportional 



to t^^'^ and the other to t ^. According to Eq. ( [BIOD , we may also use the residual gauge 
freedom for Ht to add any function which does not depend on t. One may choose it as 

Ht = if'' - ) ^{x) + {r' - to 1) ^x), (3.13) 

so that Ht vanishes at the initial time t = to- The initial density field p(a;,to) is also 
expressed by the metric. If the density contrast S = {p — pb)/ pt is sufficiently small, we 
obtain from Eqs. ( ^.7[ ) and ( ^.61 ) 

5 = -V^(t^/^^ + r^<l>). (3.14) 

Therefore, the initial density can be related with the initial metric perturbation as 

p{x, to) = ptito) [l - V (to^'^ + to . (3.15) 

The above expression is used only to relate the initial density fiuctuation to the initial metric 
perturbations in linear regime. Once the initial seed is given, the later non-linear evolution 
is evaluated by the non-perturbative expression Eq. (|2.14|) . 



B. Vector perturbations 

The Einstein equations for the vector perturbations are 

V^(/..-l6.)=4@\, (3.16) 

hi + 3-hi - —^bi = 0. (3.17) 

Introducing [3i as 

hi{x) = V2/?.(x), (3.18) 

Eq. (|3.16|) is solved to give 

= -3 (t-V3 _ t-V3) v2^^(^) ^ 8 _ ^^^^^^ ^g^^g^ 
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where we again used the residual gauge freedom (cf. Eq. ( [Bll| ) in the Appendix B) to set 
hi{x,tQ) = 0. The Pi{x) is directly related to the initial value of the vorticity field. From 
Eq. ( p.lTp , we have 



u;\x, to) = I ^ g^^''V^/3,,fc, (3.20) 
2\/-^(a:;,to) 



where e^^^ is the 3-dimensional Levi-Civita symbol with e^^^ = 1. Particularly when the 
deviation from the background is sufficiently small, we obtain the first-order expression 

u\x,to)^—^^e'^'V'P,,u. (3.21) 

C. Tensor perturbations 

The equation for the tensor perturbations is 

Hi, + 3-ify - \v^Hi, = 0. (3.22) 

This is a homogeneous wave equation in the expanding universe. The solutions are well 



known and we will not discuss the detail here. See, e.g., . 

Before closing this section, we emphasize the following point: As for the metric, our 
result is identical to that of the linear perturbation theory (e.g. |11|). However, there is an 
important difference. In Bardeen's paper, it is essential to linearize the density contrast. On 
the other hand, our Lagrangian approach does not rely on the assumption that the density 
contrast should be small. It is actually an important advantage of the Lagrangian approach 
that it uses (or extrapolates) the well-known solutions of the linear theory to express the 
non-linear density contrast. 

IV. SUMMARY AND DISCUSSION 

Motivated by the fact that the Lagrangian condition enables us to obtain simply the 
integrals of the density and the vorticity along the fluid flow, we have developed the La- 
grangian perturbation theory in the general relativistic cosmology, fully using this condition. 
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The main advantage of the present Lagrangian theory is to be amenable to the vorticity, 
while previous works are limited within the irrotational fluid. In this approach, only the 
metric is expanded and its behaviors are determined perturbatively through the evolution 
equations, while the density and the vorticity are calculated non-perturbatively from the 
integrals along the flow. Seeds of the density and vorticity fluctuations can be related with 
the metric perturbation at an initial time, only when the constraint equation is solved as the 
Poisson-like equation. Hence, this simple approach is the first that can describe the dynam- 
ics of the dust fluid with the rotational motion as well as the non-linear density fluctuations, 
up to the caustic formation. 

As an illustration, the first order Lagrangian perturbation has been solved explicitly. It 
is natural that at the first order of the metric perturbation, our result agrees with that of 
the usual linear theory when we assume the small density contrast. Contrary to the 
standard linear perturbation theory, it is explicitly shown in our approach that the vorticity 
is coupled to the density contrast and is amplified in high density (i.e., collapsing) region. 

In this paper, we have presented a relativistic framework to describe non-linear evolution 
of the density fluctuation which is not restricted to the irrotational case. Once the framework 
is given, it should be straightforward to apply it to higher-order computations. Actually, 



an extension to the second order is being investigated by Morita pj]. More comprehensive 



study of non-linear couplings will be of the subject of future investigation. 
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APPENDIX A: VORTICITY IN THE NEWTONIAN COSMOLOGY 



In this Appendix, we summarize the treatment of vorticity in the Newtonian cosmology. 

Using expanding coordinates r = a{t)x, the continuity equation and the equation of motion 
are 

dp a 1„ , ^ , , , 

-£ + 3-p+-V-{pv) = 0, Al 
ot a a 

dv a 1 , 1 „ , , . ^, 

— + -v + -{v- V)v = — V0, (A2 
at a a a 

where p is the density of dust matter, v is the pecuhar velocity defined by 

dr a dx ... 
V = — r = a—, (A3j 

dt a dt 

and is the gravitational potential which satisfies the Poisson equation 

VV = 4vrG'a^(p-p6). (A4) 

Let us introduce the Lagrangian time derivative 

d d I ^ 

dt ot a 

and the Lagrangian coordinates q = x{tQ). In the Lagrangian description, the continuity 
equation ( |A1|) is 



t + 3^p + ^V.. = 0, (A6) 
dt a a 



which is solved to give 



We now discuss the vortical part of the fluid flow. The vorticity is defined by 

u} = -V XV. (A8) 
a 

From Eq. ([A^), we obtain the propagation equation for the vorticity: 
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duj a uj „ 1 , , . , 

— + 2-U + —V ■ V = -ioj ■V)v. A9 

at a a a 

Using Eq. ([A^), we accordingly have the diffusion equation of Beltrami: 

at \ap J \ap J 

It is instructive to define the vector components of the vorticity in the Lagrangian frame as 
follows: 

. It^. (All) 

Then, Beltrami's equation ( |A1CI| ) simply reads 

0, (A12) 



d j uj 

dt \ap 



which is manifestly solved to give 



ap ap 



to 



(A13) 



or 



aopo dqi 

The result Eq. ( [A14| ) is known as Cauchy's integral 



u;\q,t) = ^^u^{q,to). (A14) 



APPENDIX B: RESIDUAL GAUGE FREEDOM IN THE LAGRANGIAN 

CONDITION 

The general gauge transformation to the first order is induced by the infinitesimal coor- 
dinate transformation 

= x^ + i^". (Bl) 
The changes due to the gauge transformation are 



11 



h9iiv — —gtiv,aC' — 9lJ.a^'^,v — gua^'^,n, (B2) 

5^M^ = e%w^-<.r. (B3) 

For the spatially flat background, are decomposed into each mode, 

e^{T,S'^Lj+e), (B4) 

where the vector mode quantity t satisfies — 0. In order to maintain the Lagrangian 
condition u^^ — {1, 0, 0, 0), we have S^ui^ — 0, which leads to ^'*^o — 0. Therefore, 

T^T{x), Li^Liix), h^h{x). (B5) 

The changes of the metric tensor due to the residual gauge freedom are 

= Ti, (B6) 

S^Qij = -a^ (^2^T5,,- + 2Lij + l^j + . (B7) 

Hence we obtain 

B^B + T{x), (B8) 

Hl^Hl--T{x), (B9) 

Ht^Ht-L{x), (BIO) 

hi^hi-ii{x). (Bll) 

The vector mode quantity 6j and the tensor mode quantity are gauge- invariant. 
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